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ABSTRACT 

In  one  method  of  solution  to  the  problem  of  making  up  a  schedule 

of  work  in  a  "job-shop"  type  industry  there  arises  a  symmetric  matrix 

the  conflicts  among  the  jobs  to  be  done.  This  thesis  is  primarily 
concerned  with  the  transformation  of  this  "conflict"  matrix  to  another 

matrix a  diagonalized  symmetric  matrix.  That  is  a  matrix  having  its 

nonzero  elements  bunched  or  clustered  near  the  principal  diagonal. 
The  transformation  is  considered  as  the  product  of  simple  transformations 
each  of  which  interchange  rows  or  columns  of  the  undiagonalized  or 
partially  diagonalized  matrix. 

A  numerical  technique  is  devised  for  the  systemmatic  reduction 
of  an  undiagonalized  symmetric  matrix  to  successively  better  diagonalized 
matrices .  The  technique  is  programmed  for  use  on  the  NCR  102A  Digital 
Computer. 
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CHAPTER  I 
INTRODUCTION  TO  THE  PROBLEM  UNDER  CONSIDERATION 

1 ,  Introduction 

In  a  "job-shop"  type  industry  there  are  many  jobs job  orders 

to  be  done,  each  requiring  the  use  of  certain  of  the  varied  facilities, 
special  tools,  and  specially  trained  technicians  available  in  the 
industry.  The  question  of  how  to  make  up  the  time  schedule  of  work 
to  be  done  arises.  That  is  which  jobs  should  be  time  scheduled  first 
considering  the  complexity  and  use  of  varied  facilities,  etc.,  in  the 
industry.  Then  having  time  scheduled  the  first  one,  which  job  should 
be  time  scheduled  next,  and  so  forth. 

The  mathematical  problem  to  be  considered  in  this  thesis  will  assist 
in  the  determination  of  an  orderly  procedure  for  time  scheduling  the  jobs 
to  be  done  and  will  provide  the  order  in  which  jobs  should  be  considered 
for  time  scheduling. 

2.  The  Mathematical  Problem 

Given  a  symmetric  matrix  A,  consisting  of  zero  and  nonzero  elements; 
find  a  symmetric  matrix  C  whose  elements  satisfy  the  f -allowing: 

1)  The  elements  of  a  given  column  of  A  become  the  elements  of 
some  column  of  C . 

ii)  The  elements  of  a  given  row  of  A  become  the  elements  of  some 
rows  of  C. 

iii)  The  nonzero  elements  of  the  matrix  C  are  located  along  the 
principal  diagonal  of  the  matrix.  That  is  to  say  there  is  a  clustering 
of  elements  near  the  principal  diagonal. 


iv)  In  addition  be  able  to  enumerate  the  correspondence  of  columns 
of  the  matrix  A  with  columns  of  the  matrix  C. 

The  conditions  i  and  ii  together  are  equivalent  to  the 
following  I 

i' )  Columns  may  be  rearranged,  maintaining  the  order  of  elements. 

ii')  Rows  may  be  rearranged,  maintaining  the  order  of  elements. 
For  consider  any  element  of  the  matrix  say  (i,  j).  The  column  i  is 
shifted  to  say  the  h-th  column  not  necessarily  retaining  the  order  of 
elements .  Also  row  j  is  shifted  to  say  the  k-th  row  not  necessarily 
retaining  the  order  of  elements.  The  element  that  was  in  position 
(i,  j)  must  now  be  in  position  (h,  k)0  However,  the  same  result  is 
obtained  by  shifting  the  i-th  column  to  the  h-th  column  maintaining  the 
order  of  elements ,  and  then  shifting  the  j-th  row  to  the  k-th  row  main- 
taining the  order  of  elements.  This  must  be  true  for  all  (i,j)#  Thus 
i  and  ii  imply  i'  and  ii1 .  That  1'  and  ii'  imply  i  and  ii  is  obvious. 

An  assumption  is  made  that  the  number  of  nonzero  elements  in  any 
row  or  column  is  small  compared  to  the  order  of  the  matrix,  say  one  tenth 
the  elements  are  nonzero.  If  this  is  not  the  case  the  idea  of  diagonal - 
ization"  or  clustering  the  elements  along  the  principal  diagonal  has  little 
significance . 

The  transformation  or  sequence  of  transformation  used  to  find  the 
matrix  C  should  be  applicable  to  matrices  of  large  order  say  62  x  62  or 
even  as  large  as  500  x  500.  Since  we  are  considering  large  order  matrices 
it  seems  logical  that  the  transformation  or  sequences  of  transformation 
should  be  susceptible  to  programming  on  a  computer,  say  the  NCR  102A 
Digital  Computer. 


3.  Notation. 
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Figure  1    General  Form  of  Matrix 


A  symmetric  matrix  A  of  order  n,  an  n  x  n  array  of  zero  and  non- 
zero elements .  The  exact  value  of  the  nonzero  elements  is  not  important — 
the  reasons  for  this  are  explained  in  Chapter  IV. 

All  nonzero  elements  P..  are  located  by  coordinates  (C. ,R.) 

ij  '  i  j 

indicating  its  position  in  the  matrix  array  is  in  the  C,th  column 
and  the  R.th  row.  All  elements  not  designated  and  enumerated  are 
assumed  zero. 

The  columns  are  designated  C,  ,C_, ,C  numbered  in  the  usual 

fashion  from  left  to  right. 


The  rows  are  designated  R,  ,R  , R  numbered  from  the  bottom 

upward.  This  unusual  numbering  is  associated  with  the  numerical  technique 
used  to  "diagonalize"  the  matrix. 

The  sum  of  the  row  and  column  positions  is  designated  S. . ,  where 
S . .  =  C.+  R. .  Note  that  all  the  elements  of  the  principal  diagonal 
have  the  same  value,  namely,  n  +  1. 

The  distance,  D,  of  an  element  from  the  principal  diagonal  is 
S. .  -(n  +  1)  =  D.  This  is  an  integer  distance  along  the  ith  column  or 
the  jth  row. 
4-.  The  Measure  of  Diagonal  izat ion  of  a  Matrix 

Some  criteria  must  be  set  up  to  measure  the  clustering  or  bunching 
of  the  nonzero  elements  along  the  principal  diagonal  to  distinguish 
between  the  matrices .  That  is  to  determine  which  of  two  matrices  is 
more  diagonal ized.  This  measure  is  relative  and  serves  only  as  an 
index. 

A  simple  criterion  would  be  to  minimize  the  distance  of  the  element 
farthest  from  the  principal  diagonal. 


Min 


Max 


(n  +  1)  -Si.   {     Min  Max  D 


This  will  produce  a  certain  amount  of  diagonalization  ,  However,  a 
simple  example  will  show  that  a  more  powerful  criterion  is  needed. 
Figures  2a  and  2b  both  have  the  same  Min  Max  D,  namely  2,  however, 
Figure  2b  is  more  compact  along  the  principal  diagonal. 
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Figure  2  Examples  of  Min  Max  D  Criteria 


Consider  the  criterion  that  the  sum  of  the  distance  from  each 
element  to  the  principal  diagona]  squared  be  minimized. 
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Each  element  far  from  the  principal  diagonal  will  contribute  greatly 
to  the  summation,  where  is  elements  next  to  the  principal  diagonal  will 
have  negligible  effect  on  the  summation.  Thus  the  greater  the  cluster- 
ing along  the  principal  diagonal,  the  smaller  will  be  the  summation. 
The  previously  considered  Min  Max  D  criterion  will  tend  to  be  satisfied 
also .  For  consider  Figure  3  showing  two  columns  (all  other  columns  are 
assumed  the  same  in  both  cases  and,  therefore,  contribute  the  same  to 
the  summation  in  both  cases). 
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Figure  3 


(b) 
Example  of  /  D  for  Row  Interchange 


In  the  summation  note  that  only  the  contribution  of  the  outer  element 
in  each  column  is  affected  by  the  interchange  of  the  columns  under  con- 
sideration in  Figure  3.  In  (a)  the  outer  elements  contribute  to  the 
summation  (k+1)  +  (k-l)     or  2k  -t-2.  In  (b)  these  same  elements  con- 
tribute  2k  .  Thus  (b)  is  better  diagonalized  in  the  sense  of  both 
Min  Max  D  and  Min  >  D2. 


For  the  remainder  of  the  discussion  and  in  testing  the  merits  of 


a  transformation  under  consideration  we  will  use  Min  <-_D_. 


CHAPTER  II 
SOME  FEATURES  OF  MATRICES  AND  ELEMENTARY  TRANSFORMATIONS 

A  square  matrix  A  whose  elements  satisfy  the  conditions 

aiJ  =  aji 
possesses  symmetry  about  its  principal  diagonal,  and  is  called  a 
symmetric  matrix. 

If     C  =  B  A,   then  the  matrix  A  is  said  to  be  linearly  trans- 
formed into  the  matrix  C.  The  matrix  B,  which  affects  the  transformation 
of  A,  may  sometimes  be  advantageously  thought  of  as  the  multiple  product 
of  very  simple  matrices. 

B  =  B  B  .  B.  B, 

n  n-1       2  1 

The  transformation  of  A  is  thus  accomplished  through  a  succession  of 
simpler  transformations,  the  first  of  which  is  B.  ,  the  second  B_,  and 
so  forth. 

The  matrix  can  be  transformed  through  postmultipllcation  as  well  as 
by  premultiplication,  i.e.,   C  =  A  B' .   In  which  case  the  matrix 
affecting  the  transformation  is  decomposed  into  components 


2 Bn 


B'  =  B1  B B. 


Where  E,  indicates  the  first  simple  transformation,  B?  the  second,  and 
so  forth. 

It  is  of  interest  to  determine  the  simplest  fundamental  forms  into 
which  can  arbitrary  nonsingular  matrix  B  may  be  decomposed,  and  to 
interpret  the  individual  transformations  which  they  effect  upon  the  form 


of  the  matrix  A.  These  are  called  elementary  transformations  and  the 
matrices  which  produce  them  are  called  elementary  transformation  matrices. 
There  are  three  fundamental  types  of  elementary  transformations  and 
correspondingly  three  fundamental  types  of  transformation  matrices . 
The  first  type  amounts  to  an  interchange  of  two  rows  or  columns  of  the 
matrix  A.  The  second  type  is  the  addition  of  the  elements  of  a  row  or 
column  to  the  corresronding  elements  of  another  row  or  column.  The  third 
type  is  the  multiplication  of  any  row  or  column  of  A  by  an  arbitrary 
nonzero  factor.  For  our  purpose  it  is  necessary  to  develop  only  type 
one  to  a  greater  extent. 

The  transformation  may  be  affected  by  multiplying  A  by  a  transforma- 
tion matrix.   If  the  desired  transformation  is  intended  to  affeqt  the 
ROWS,  A  is  PREmultiplied  by  the  transformation  matrix;  if  the  COLUTfJS 
are  to  be  affected,  A  is  POSTmultiplied  by  the  transformation  matrix. 

The  elementary  transformation  matrix  is  formed  from  the  unit  matrix 
by  performing  on  it  the  same  elementary  transformation  which  the  desired 
elementary  transformation  matrix  is  intended  to  effect  in  the  matrix  A  by 
means  of  (pre-)postmultipJLication.  Thus  I  (the  unit  matrix)  with  its 
pth  and  qth  rows  (columns'  interch  mged  yields  a  transformation  matrix 
which  by  means  of  pre-  (post-Nmultiplication  interchanges  the  pth  and  qth 
rows  (columns^  of  the  matrix  A.  Apparently  the  transformation  matrix 
has  two  forms  depending  on  whether  it  is  intended  to  effect  a  transforma- 
tion of  the  rows  or  columns  of  A,  however,  it  is  to  be  noted  that  one 
form  is  the  transpose  of  the  other,  and  further  for  type  one  the  two  forms 
are  identical.  */e  will  denote  this  elementary  transformation  matrix  as 

T    where  the  subscripts  indicate  the  rows  (columns)  being  interchanged, 
p-q 


10  0  0' 

T 
2-3 

_     0  0  1  0! 

0  1  0  O1 

=       T 

2-3 

Rows 

[0  0  0  ll 

Columns 

As  an  example ,  for  a  fourth  order  matrix  to  interchange  the  second  and 
third  rows  (columns  \ 


T      =    T 
23        32 


Consider  a  symmetric  matrix  A,  pre-  and  postmultiplied  by  the 
elementary  matrix  T.  ..  The  resulting  matrix 

A'  =  W  Tn 

is  a  symmetric  matrix  with  both  the  ith  and  jth  rows  and  columns  inter- 
changed .  Continuing  this  process  we  may  by  the  elementary  matrix 

T,   produce  the  matrix 
K-m  r 

k' '   =  W  Tk-M  =  Tk-m  Ti-J  A  Tl-J  Tk-m 

Thus  by  a  multiple  product  of  elementary  transformation  matrices  pre- 
and  postmultiplying  the  matrix  A,  the  matrix  A  is  linearly  transformed 
into  a  symmetric  matrix  C. 

The  matrix  C  which  may  be  obtained  from  the  matrix  A  bv  means  of  a 
finite  number  of  elementary  transformations  is  said  to  be  equivalent  to  A, 
The  equivalence  of  matrices  is  a  mutual  relationship  since  the  trans- 
formations are  nonsingular  and  hence  reversible .  Thus  if  C  can  be 
obtained  from  A  by  a  succession  of  elementary  transformations,  it 
follows  that  A  can  be  obtained  bv  means  of  elementary  transformations 
from  C .  The  state  of  equivalence  of  the  two  matrices  A  and  C  may  be 
stated 


P  A  Q  =  C        where  P  and  Q  are  nonsingular  matrices . 

'  ow  if  A  is  a  symmetric  matrix  arrl  Q  =  Q,  Q_ Q  ,  where  each 

3,  is  an  elementary  transformation  matrix  T,  if  further 

p  =  aT=  (;1q2....^)t  =  qd *& 

then  C  is  also  a  symmetric  matrix.  Since 

T 
C  =  Q  A  Q 

The  transpose  of  C  is  then 

C  =  Q  A1  (Q1)  =  Q1  A  Q. 

A  matrix  which  is  its  own  transpose  is  symmetric .  So  C  is  symmetric . 


The  preceding  discussion  has  covered  the  elementary  transformation 
matrices  in  general .  Now  to  apply  these  concepts  in  the  notation  to  be 
used  in  the  numerical  technique  to  be  developed  in  Chapter  III. 

In  order  to  have  a  symmetric  matrix  after  premultiplication  by 
an  elementary  transformation  matrix  T.  .  which  interchanged  the  ith  and 
jth  columns  we  must  multiply  by  an  appropriate  row  interchanging  elemen- 
tary matrix  transformation.  Since  the  rows  are  designated  from  the 
bottom  up,  we  shall  designate  the  appropriate  postmultiplying  elementary 

matrix  as  T — - ,     Where 
i  J 

1  ♦  i  =  (n  +  1) 

j  *   J  =  (n  ♦  V 

10 


Examination  of  Figure  1  will  show  that  this  indeed  is  the  appropriate 

row  interchanging  matrix. 

Thus  AJ  =  (T±.   A)  T— . 

As  an  example  consider  the  matrix  A  interchanging  first  the  2nd  and 

5th  columns,  then  the  corresponding  rows. 


A  = 


1  0 

0  1 

1  1 

0  0 

1  1 


A'  = 


1110  0 
0  110  1 
10  111 
0  0  110 
110  0  1 


*1  = 


1110  0 
110  0  1 
10  111 
0  0  110 
0  110  1 


then  i 
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A'    =  T25  A 

6-2=4- 

6-5=1 


A1  -  A'   T^       «     *u 


25 


=  A«   T, 
A  )  T, 


ka 


Similarly  it  can  be  seen  that  if  the  postmultiplication  elementary 
transformation  matrix  is  performed  first  then  followed  by  the  appropriate 
premultiplication  elementary  transformation  matrix  we  have 


*2  =  TJJ  (A  T   ) 


As  an  example  consider  the  matrix  A  interchanging  first  the  2nd  and 
5th  rows ,  then  the  corresponding  columns . 


10  10  1 

0  110  1 

=11110 

0  0  1  1  oj 

1  1  0  0  lj 


A'    = 


0  0  110 

0  110  1 

11110 

10  10  1 

110  0  1 

A2  = 


10  10  0 
0  110  1 
11110 
0  0  111 
0  10  11 


11 


A-  =  A  T25 

A2  =  TU  A'  =  TU  V 

Note  that  T   applied  as  a  premultinlication  matrix  and  as  a  post- 
multiplication  matrix  produces  different  transformed  matrices.  That  is, 

TU  (  AT25  >  ■*a  +  *l"(T254)t41 

Since  we  are  discussing  transformations  of  symmetric  matrices  only, 
necessarily  these  elementary  matrix  transformations  will  occur  in  groups 
of  two.  So  for  convenience  define  a  pair  of  elementary  matrix  transforma- 
tions as  the  transformations  T. .  applied  to  the  rows (columns)  and  the 
corresponding  T7T  applied  to  the  columns  (rows\ 

The  given  matrix  A  has  a  certain  <>_  D  ,  say  (  >__    D  )_,  which  in  the 
previous  example  was  £2.  If  we  apply  the  transformations 


*i  =  (Tu  A'  TU 


sr-     2 
the  transformed  matrix  will  have  a  new    ^  D',  which  for  the  previous 

example  is   JP .     If  on  the  other  hand  we  apply  the  transformations 

A^Tjj   (A  Ti}) 

a  different  transformed  matrix  is  obtained  with  a  certain  ^_  D  which 
for  the  previous  example  is  32.  Thus  it  can  be  seen  that  different  pairs 

of  elementary  transformation  matrices  will  produce  transformed  matrices 

it  ii 

diagonal ized  to  a  greater  or  lesser  degree.  Some  pairs  of  elementary 

matrix  transformations  may  even  produce  a  transformed  matrix  which  is 

less  diagonal! zed  in  the  sense  of  Min  ^D   . 


12 


CHAPTER  III 

PROCEDURE 

1 .  General 

Basically  the  procedure  is  to  apply  a  pair  of  elementary  trans- 
formation matrices  to  the  undiagonalized  or  partially  diagonalized  matrix 
and  test  the  diagonalization  in  the  sense  of  minimizing  ^_  D  . 

a.  If  the  transformed  matrix  is  better  diagonalized,  then  use  this 
matrix  as  the  matrix  to  be  diagonalized  and  repeat  the  process  applying 
elementary  transformation  matrices  to  this  partially  diagonalized  matrix. 

b.  If  the  transformed  matrix  is  NOT  better  diagonalized,  apply  a 
different  pair  of  elementary  transformation  matrices  to  the  undiagonalized 
matrix  and  test  for  ^__  D  again. 

We  can  choose  the  pair  of  elementary  transformation  matrices  either 
by  some  random  process  or  systematically.  This  thesis  is  concerned  with 
a  systematic  procedure  for  improving  the  diagonalization  of  a  given 
matrix. 

2.  A  systematic  Procedure 

An  element,  let  us  designate  it  P,,,  which  lies  farthest  from  the 
principal  diagonal  individually  contributes  the  most  to  2L.D  of  any 
element  P..  in  the  matrix.  Since  by  symmetry  there  are  two  such  elements 
let  us  consider  the  one  with  the  largest  S,  that  is  STJ.  It  appears  that 
a  pair  of  elementary  matrix  transformations  which  involve  this  element 
will  be  a  good  transformation  to  hypothesize  making.  Furthermore,  to  move 
this  element  to  a  position  one  unit  from  the  principal  diagonal  will 
reduce  the  contribution  to  ^D  of  this  element  to  a  minimum.  However, 
other  elements  of  the  rows  and  columns  interchanged  will  change  distance 

13 


from  the  principal  diagonal  so  that  one  must  test  this  hypothesized  trans- 
formation to  see  if  it  improves  the  diagonalization.  If  it  does  we  make 
the  hypothesized  transformation  and  repeat  the  process  again  choosing 
the  new  P--  with  a  maximum  S,,. 

If  the  previously  hypothesized  transformation  does  not  improve  the 
diagonalization  we  hypothesize  another  transformation.  If  the  previous 
elementary  transformation  pair  T.  .  was  a  tremultiplication  first,  i.e., 
row  interchange  first,  we  now  hypothesize  T,.,  making  a  column  interchange 
first  such  that  P,.  is  moved  to  one  unit  from  the  principal  diagonal.  If 
this  results  in  a  better  diagonalized  matrix  we  make  this  hypothesized 
transformation  and  repeat  the  process  again,  choosing  the  new  PTJ  with  a 
maximum  STJ. 

If  the  hypothesized  transformation  moving  P,,  to  a  position  one 
unit  from  the  principal  diagonal  does  not  improve  the  diagonalization  of 
the  matrix,  hypothesize  moving  P_-  to  two  units  from  the  principal 
diagonal.  If  neither  of  these  produce  a  better  diagonalized  matrix, 
hypothesize  moving  P....  to  three  units  from  the  principal  diagonal  and 
continue  until  the  interchange  is  with  the  row  (column)  adjacent  to  the 
element  PTT. 

If  all  these  hypothesized  transformations  with  this  element  P,,  fail 
to  produce  a  better  diagonalized  matrix  disregard  this  element  and  choose 
from  among  the  remaining  elements  the  element  with  the  largest  S . . .  Repeat 
the  process  of  hypothesizing  transformations  with  this  element  until  one 
is  found  which  improves  the  diagonalization.  If  none  is  found  among  all  the 
elements  with  S  .  greater  than  n  +  1,  i.e.,  we  have  looked  at  all  possible 
moves  of  all  elements  on  one  half  of  the  symmetric  matrix,  then  the  matrix 
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which  we  are  attempting  to  diagonalize  is  itself  the  best  possible 
diagonalization  in  the  sense  of  minimizing   <£_D;->. 

By  hypothesizing  moving  this  "wildest"  element  to  one  unit  from 
the  principal  diagonal  first,  we  are  maximizing  the  possible  improvement 
in  on--  stage  of  diagonalization.  The  actual  improvement  achieved  in  one 
stage  will  usually  be  somewhat  less  however.  The  distance  between  the 
rows  (columns)  being  interchanged  must  be  variable  since  it  can  be  shown 
that  a  system  which  hypothesizes  making  the  interchange  always  with  an 
adjacent  row  (column)  will  not  achieve  maximum  diagonalization. 
Figure  4a  is  the  best  diagonalization  possible  using  a  hypothesis  rule 
of  always  making  the  interchanges  with  the  adjacent  column  (row). 
Figure  4-b  is  the  best  possible  diagonalization  using  the  systematic 
procedure  outlined  above  on  the  same  matrix. 
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Figure  4-    Example  of  Systemmatic  Procedure  For 
Better  Diagonalization 


3.  Systematic  Procedure  Applied  to  the  NCR  102A  Computer 

All  column  and  row  numbering  is  done  in  the  octal  number  system 
for  the  simple  reason  that  the  NCR  102A  computes  in  binary  number 
system  using  data  inserted  in  the  octal  number  system.  This 
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OVERALL  FIOW  CHART-DIAGONALIZATION  PROCESS 
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Figure  5 
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information  could  be  inserted  as  decimal  numbers.  However,  before 

the  computer  can  start  to  compute  a  subroutine  would  have  to  be 

employed  to  convert  all  data  to  octal. 

The  form  and  amount  of  information  about  a  particular  nonzero 
element  of  the  matrix  A  which  will  be  placed  in  the  computer  is  dictated 
in  part  by  the  computer  and  in  part  a  desire  by  the  writer  to  contain  in  one 
computer  storage  word  all  information  relative  to  this  particular  element. 
Since  a  computer  storage  word  is  composed  of  12  octal  digits  it  is  conven- 
ient to  store  the  information  about  a  particular  P. .  as  follows 2 
i)  for  matrices  of  order  63  x  63  or  less 


C©     ftjo         ^>cj      Cj      Kj 


Since  numbers  as  large  as  63  can  be  represented  in  2  octal  numbers  (63 

decimal  is  77  octal)  this  representation  is  possible.  Where 

C.    Column  designation  of  the  nonzero  position  P. .  in  the  matrix  before 

any  diagonalization  has  been  effected. 
R    How  designation  of  the  nonzero  position  P. .  in  the  matrix  before 

any  diagonalization  has  been  effected . 
C.    Current  column  designation  of  this  P. .  is  the  matrix  at  any  stage 

of  diagonalization. 
R.    Current  row  designation  of  this  P. ,  in  the  matrix  at  any  stage 

of  diagonalization. 
S.  .   The  current  C .  ♦  R . . 

It  is  to  be  noted  that  the  C. ,  R  ,  and  S, .  may  change  from  time  to  time 

J  v 

during  the  diagonalization  process.  However,  the  C.  and  R.  remain 
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unchanged  and  thus  a  one  to  one  correspondence  between  the  rows  (columns) 
of  the  original  matrix  and  the  current  partially  or  completely  diagonal- 
ized  matrix  is  obtainable  by  comparing  the  R.  and  R,   (C.  and  C.) . 
The  detailed  numerical  techniques  described  later  are  based  on 
stored  word  structures  for  order  63  x  63  or  less . 

ii)  for  matrices  of  order  511  x  511  or  less 


O  O 


■   .v  ■  - 


y  , ,  i l 


Ceo     ^0    Cc      R^ 


Numbers  as  large  as  511  can  be  represented  in  3  octal  numbers  (511  decimal) 
is  777  octal \  However,  certain  information  must  be  sacrificed  in  order 
to  contain  the  essential  information  in  one  computer  storage  word  length. 
For  these  larger  matrices  the  value  of  S. .  must  be  formed  wherever  needed 
rather  than  being  immediately  available.  The  symbols  have  the  same 
meaning  as  above. 

Determination  of  the  Largest  S... 

In  order  to  determine  the  "wildest"  element  or  more  precisely  the 
element  with  the  largest  S.  .,  it  is  necessary  to  order  the  elements  P.. 
or  at  least  order  some  of  them,  say  the  largest  6.  The  reason  that  more 
than  one  largest  S. .  must  be  obtained  is  that  it  may  happen  that  none  of 
the  hypothesized  transformations  using  the  P..  element  with  the  largest 
S,j  will  improve  the  diagonalization.  However,  it  is  unlikely  that 
among  the  elements  with  the  6  largest  S.  .!s  at  least  one  will  not  have  a 
transformation  which  will  improve  the  diagonalization.  If  one  of  them 
does  improve  the  diagonalization  the  resulting  Pj*'*  will  be  reordered 
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Figure  6 
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and  the  6  largest  S. .'s  of  the  new  matrix  obtained.  It  is  to  be  noted 
that  the  technique  of  ordering  tends  to  order  not  only  the  6  largest  but 
to  relatively  order  all  the  elements.  Thus  after  several  partially 
diagonalized  matrices  are  obtained,  all  the  P. .'s  will  be  relatively  well 
ordered . 

The  technique  used  is  simply  to  compare  the  S.  .  of  the  K  and  K  +  1 
storage  cells,  putting  the  P. .  with  the  smaller  S. .  in  the  K  storage  cell. 
Then  compare  the  S. .  of  the  K  +  1  and  K  +  2  storage  cells,  putting  in 
K  +  1  the  P. .  with  the  smaller  S. ..  Continuing  this  process  until,  the 
N  (last)  storage  cell  is  reached,  then  put  both  P*j's  away.  Decrease 
the  number  of  storage  cells  under  consideration  by  one,  that  is,  for  the 
second  time  through  this  ordering  process  only  the  first  N-l  storage  cells 
woul :  have  their  S   's  compared.  To  determine  if  the  preassigned  number 
of  largest  S  's  have  been  obtained,  test  if  N-6  is  the  last  storage  cell 
to  be  considered  for  the  next  time  through  the  ordering  process.  If  not 
repeat  the  process  starting  with  the  first  two  storage  cells ;  if  so  go 
on  to  the  next  segment  of  the  program. 


Inking  a  Hypothesized  Transformation 

For  the  first  hypothesized  transformation  we  consider  interchanging 
the  elements  of  two  rows  first.  The  P_T  stored  in  the  N  storage  cell  will 
have  the  largest  STT.   It  is  easily  determined  which  row  this  PTT  is  in, 
namely  the  Jth  row,  and  the  distance  from  the  principal  diagonal,  namely  D, 
Then  J-D  is  the  row  which  has  the  Ith  column  on  the  principal  diagonal. 
Let       K  =  J  -  D  +  1 

Then  the  interchange  of  the  Jth  and  Kth  rows  will  place  PT,  next  to  the 
principal  diagonal .  Now  to  maintain  symmetry  two  columns  must  also  be 
interchanged .  It  is  easily  seen  that  the  columns 

20 


FJ^W  CHART  -  MiKDJC  A  HYPOTHESIZED  TRANSFORMATION 
— _ _^ 


© 


C    .- 


Find  P   In  Xth  storage  cell 

Print  P  used  to  determine  transformation 
Form  I  =  S  -  (n  ♦  1}        


Form  in  buffer: J,  K=J-Y«-1  in  last 

2  octal  placess 
Form  in  buffer:  JHn^-.", and  J=(n+1)-J 

shifted  left  2  octal  places 
Initialise^ lookup  and  put  away  for  P 

Place  P  in  buffer 
(Form   (J-J)  and  (K-jJ 


-<K-J)>0  /*Yes  — 


'(j-j;>0  — »Yes 

No 

I  Extract"  j-K  i:, 

|  Form    J-i  and  K-i  [ 


(j-i)>0>»Yea 


I 

No 
I 
[Extract  j=3  into  P~] 


4 

No  No 

(lxtr.act.irS.ir.to  pj     (Extract  I 

porm  new  '" .  extract  into~P  "         | 

Form  D=f.-'n«-l" 
IForm  D2,  -         -;al  Zdj  | 

j  Put  away  P 

!  Alter  lookup  ->>sses  | 

J\    /'  —I«    '-'+1  >  address  of  P? 
Yes«— <JM  10> — *No 


Print 


,  _c 


-No<5 


Reload  memory  with 
partially  ordered 
P|_s  previously  saved 


*Yes 


Is  previous     X  D2 
>this    5_  D2 


Form  in  buffer:  I,   L=l-Y*l   shifted 

left  2  octal  places 
Form  in  buffer!  L=(n+ll-L,I=(n+l~-I 

in_.last_2_  octal_pl.aces 


[initialize  lookup  and  put  away  for  P~] 

ri^ce  P  in  buffer 
Farm   (I-ll,   and    (W) 


*Yes  - 


Extract  1=L  into  pi 

:Tr.:T_ 


-<(L-i  )>0) ^Yes  - 

fextrac t  1=K  into  P  j 


[Form  (I- J )  and  (L-J) 


>Yes- 


T 


>Yes 


No 
[Extract  j=K  into  p| 


|Extract  J=L  into  P  | 

?urm  new  S,  I    into  P 

Form  Da6-(n+]  ' 

gorra  D2,    idd  '  .1    2.  >- 

1 
|Pi:t   luay  P 

,'Alter  lookup  and  ;  ut   iwa\ 

,'    v^'"'"2 Is  K*^  >  address  of  P 

Yes<—    TM  l2>-*No  , 

V"'  j    Print    S  D2 


;] 


Yes*" 


WNTM  13     • 


'--If 


Tally  in  iteration 
^int  Iteration  and£D2 

j  Recycle 


<5 


■^■So 


previous  £D2 
this  JT  D2 


a 


Reload  memory  with 
partially  ordered 
PJ s  previously  saved, 


<$-L>N/fc> >Yes  — >(^T) 

4 

.0 

rzib~] 


(pondence  of  Original  "1 
lew  Columns  routine 


To  Correspondence 
an.i  new 


0 


J  =  (n  ♦  1)  -  J 

K  =  (n  ♦  1)  -  K 
are  the  appropriate  columns  to  interchange  to  maintain  the  symmetry. 
These  interchanges  are  the  equivalent  of  the  elementary  matrix  trans- 
formation 

TJK  A  T  JK. 

Forming  these  quantities  J,  K,  J,  K  where  they  are  readily  available 
proceed  to  examine  each  P. .  performing  the  following  tests : 
i)  if  j  =  J  then  make  j  =  K 
if  J  =  K  then  make  j  =  J 
otherwise  row  designator  remains  unchanged . 
ii)  if  i  =  J  then  make  i  =  K 
if  i  =  K  then  make  i  =  J 

otherwise  column  designator  remains  unchanged. 

2 

Determine  the  new  S. .  of  this  element  and  then  D.  Form  D  and  add  the 

contribution  of  this  element  to  the  partial  2-  D  already  obtained. 

When  all  the  P^'s  have  been  so  treated,  test  to  determine  if  the 

<=— -   2  >c-  2 

? .  D  for  the  hypothesized  change  is  less  than  ,>_D  for  the  matrix 

r—  2 

before  transformation.  If  indeed  this  new  ^>_  D  is  less,  we  have  a 

matrix  more  diagonalized  than  before  and  we  start  the  process  over  again, 
ordering  by  the  S. . ' s ,  etc . 

If,  however,  the  new  ^_  D  is  not  less  than  before  discard  this 
hypothesized  transformation  and  replace  the  storage  cells  with  the 
partially  ordered  set  of  P^'s  which  had  previously  been  stored  separately. 
This  time  hypothesize  a  transformation  by  interchanging  elements  of  two 
columns  first.  Taking  again  P,j  in  the  N  storage  cell  determine  which 
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column  P.-  is  in,  namely,  I,  and  the  distance  from  the  principal  diagonal. 
Then  I  -  D  is  the  column  which  has  the  Jth  row  on  the  principal  diagonal. 
Let     L=I-D  +  1 

Then  the  interchange  of  the  Ith  and  Lth  columns  will  place  PTT  next  to 
the  principal  diagonal .  Now  to  maintain  symmetry  two  rows  must  also  be 
interchanged .  It  is  easily  seen  that  the  rows 

I  =  (n  ♦  1)  -  I 

L  s  (n  ♦  1)  -  L 

are  the  appropriate  rows  to  interchange  to  maintain  the  symmetry.  These 
interchanges  are  equivalent  to  the  elementary  matrix  transformation 

Pair  TIL  A  TIL- 

Forming  these  quantities  I,  L,  It  L  where  they  are  readily  available 
proceed  to  examine  each  P..  performing  the  following  tests; 
i)  if  j  =  I   then  make  j  =  L 
if  j  =  L   then  make  j  =  I 
otherwise  row  designator  remains  unchanged, 
ii)  if  i  =  I   then  make  i  =  L 
if  i  =  L   tben  make  i  =  I 

otherwise  column  designator  remains  unchanged. 

2 

Determine  the  new  S. .  of  this  element  and  then  D.  Form  D  and  add  the 

contribution  of  this  element  to  the  partial  ^L  D  already  obtained. 

When  all  the  P. .  have  been  so  treated  test  to  determine  if  the 
ij 

2Z  D  for  the  hypothesized  change  is  less  than  Z_   D  for  the  matrix 
before  transformation.  If  indeed  this  new  ^  D  is  less,  we  have  a  matrix 
more  diagonalized  than  before  and  we  start  the  process  over  again,  ordering 

23 


by  the  S .  . ' s ,  etc . 

If,  however,  the  new  ^  D  is  not  less  than  the  ^  D  for  the  matrix 
before  transformation,  discard  the  hypothesized  transformation,  replace 
the  storage  cells  with  the  partially  ordered  set  of  P^'s  which  had 
previously  been  stored  away. 
Now  make  K'  =  K  ♦  1 
L«  =  L  +  1 
If  K'  <C  J  hypothesize  making  interchanges  as  before  using  K'  and  L' 
in  place  of  K  and  L. 

If,  however,  K'  =  J,  it  means  that  we  have  hypothesized  all  possible 
row  (column )  interchangee  which  shift  the  Jth  row  (Ith  column)  and  none 
of  them  improve  the  diagonalization  of  the  matrix.  Make  N1  =  N  -  1, 
that  is  repeat  the  hypothesizing  transformations  process  using  the  next 
wildest  element.  If  none  of  these  hypothesized  transformations  will 
improve  the  diagonalization  of  the  matrix,  make  N'  =  N  -  2,  and  so  on, 
continuing  this  process  until  all  elements  have  been  tried .  If  none  of 
the  hypothesized  transformations  improve  the  diagonalization  of  the  matrix, 

then  the  best  diagonalization  possible  using  the  criteria  of  minimize 

2 

D  has  been  attained. 


I- 
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CHAPTER  IV 
A  PHYSICAL  PROBLEM 

Consider  a  job-shop  type  industry  such  as  a  shipyard  which  may  be 
considered  typical.  The  products  may  be  considered  the  accomplishment 
of  job  orders,  each  of  which  may  require  the  use  of  special  facilities, 
special  machines  and  equipment,  and  the  special  talents  of  specialized 
technicians. 

Let  us  consider  several  typical  or  probable  job  orders: 

Job  Order  #1;  Main  Battery  Range  finder;  clean,  repaint  pris- 

matical  marks  on  lenses,  reevacuate  rangefinder  tube,  fill  evacuated 

tube  with  nitrogen,  remount  and  realign  the  rangefinder. 

Job  Order  #21  Main  Battery  Fire  Control  Equipment;  correct 

systemic  errors  in  synchro  system  between  radar  and  optical  portions 

of  fire  control  system,  also  accomplish  electronic  modification  XX. 
Job  Order  #3:  Binoculars;  collimate  and  adjust  60  pairs  of 

binoculars . 

Job  Order  #4.;  Radio  Equipment;  overhaul  model  TY  receiver. 

And  so  forth .... 

The  job  orders  are  analyzed  with  respect  to  special  machines, 
facilities,  and  specialized  technical  work  required.  For  example,  Job 
Order  #1  will  require,  1)  all  the  facilities  of  the  optical  section  for 
work  space  during  cleaning,  2)  a  specialized  technician  to  repaint  the 
prismatical  marks  on  the  lenses,  3)  special  equipment  to  evacuate  the  range 
finder  tube  afterwards,  4.)  special  equipment  to  refill  the  range  finder 
tube  with  nitrogen,  5)  special  alignment  equipment  and  technical  assistance 
during  remounting. 
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Job  order  #1 
Job  order  #2 
Job  order  §3 
Job  order  #4. 
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Job  order  #2,  while  primarily  the  concern  of  electronics  technicians 
and  requiring  electronics  equipment  also  will  require  a  synchro  expert 
and  technicians  from  the  optical  shop. 

Job  order  #3  requires  the  shop  facilities  of  the  optical  shop.  Thus 
each  job  order  can  be  analyzed  and  a  chart  prepared  as  indicated  in 
Figure  P,  where  there  are  m  job  orders  and  n  columnar  headings. 

Figure  8  in  effect  represents  a  m  x  n  matrix  of  zero  and  nonzero 
elements.  Each  x  in  the  figure  represents  a  nonzero  element.  All  the 
nonzero  elements  of  a  given  row  represent  the  columns  (men,  machines  and 
facilities)  which  are  affected  by  a  particular  job  order.  Post  multipli- 
cation of  this  m  x  n  matrix  M  with  an  n  x  m  matrix,  say  w  -  the  matrix 
transpose-  will  produce  an  m  xm  matrix  A.  The  matrix  multiplication  is 
carried  out  in  the  usual  way  using  the  conventions 

(nonzero  element)  x  (nonzero  element)  =  (nonzero  element) 

(nonzero  element)  x  (  zero  element)  =  (  zero  element). 

This  matrix  M  NT  is  symmetric  with  the  rows  and  columns  representing 
job  order  numbers.  Thus  a  nonzero  element  in  the  matrix  position  (i,j) 
indicates  that  job  order  i  and  job  order  j  have  some  conflict  between  them. 
This  conflict  may  be  that  both  job  orders  require  the  use  of  the  space 
facilities  of  the  optical  shop  or  the  use  of  the  special  nitrogen  filling 
equipment.  This  matrix  product  MM  is  thus  the  "Job  order  Conflicts  matrix. 

This  symmetrica"'  conflict  matrix  is  the  matrix  which  is  diagonal ized 
elsewhere  in  this  thesis .  The  diagonalization  provides  an  ordering  of  the 
job  orders  such  that  all  job  orders  which  conflict  with  a  given  job  order 
are  listed  in  the  ordering  near  that  given  job  order.  Thus  a  job  order 
which  conflicts  with  no  other  job  orders  will  appear  at  one  end  of  the  ordering, 
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Note  that  the  ordering  is  not  unique  for  the  ordering  may  be  enumerated 
in  reverse  order  for  the  same  diagonalized  matrix. 

Let  us  consider  again  the  original  job  order  vs  men,  machines  and 
special  facilities  matrix  M.  Since  there  is  nothing  special  about 
the  order  that  the  job  orders  are  listed,  let  us  rearrange  the  rows 
so  that  the  job  orders  which  conflict  with  each  other  are  near  each 
other.  That  is  precisely  the  ordering  of  the  job  orders  given  by  the 
diagonalization  of  the  "Job  Order  conflict"  matrix.  When  the  matrix 
with  the  rows  ordered  in  this  fashion  is  scrutinized,  it  will  be  noted 
that  there  is  a  natural  rearrangement  of  the  columns  which  will  yield 

a  matrix  compact  along  a  line for  square  matrix  this  line  is  the 

principal  diagonal  and  for  the  rectangular  matrix  this  line  runs  from 
corner  to  corner.  The  results  will  be  a  matrix  such  as  Figure  9. 

By  some  criteria  such  as  the  type  commanders  decision  that  job  A 
is  most  important,  the  first  job  order  is  selected  for  scheduling.  Other 
criteria  for  finding  the  starting  place  are  to  use  the  job  order  row  with 
the  longest  length  or  the  middle  row.  Having  selected  the  starting  row 
by  some  criteria  the  complete  time  scheduling  of  job  order  A  limits 
certain  aspects  of  job  orders  B  and  C  (see  Figure  9) .  The  time  scheduling 
of  B  and  C  in  turn  limit  D  and  E.  This  continues  until  a  point  is  reached 
where  F  and  G  are  the  job  orders  being  time  scheduled.  They  can  be 
scheduled  independently  of  each  other  limited  only  by  the  job  orders 
already  scheduled.  This  is  because  they  use  none  of  the  same  men, 
machines  or  facilities.  At  this  point  the  remainder  of  the  upper  half 
can  be  scheduled  without  concern  for  the  lower  half  causing  interference. 
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APPENDIX 

i 

The  complete  program  is  included  here  for  reference  purposes .  The 
reader  is  also  referred  to  the  Flow  Charts  for  individual  segments  of 
the  program. 

The  program  occupies  cells  1200  to  174-7  in  the  main  memory.  The 
program  may  be  operated  either  minimum  access  or  sequential  numbering 
drum.  As  written  the  program  requires  also  the  use  of  a  magnetic  tape 
having  blocks  500  to  620  available  in  both  mode  1  and  mode  2.  The  nonzero 
elements  are  stored  starting  in  cell  0000  up  to  1077.  Channel  11  of  the 
memory  is  reserved  for  temporary  storage  during  the  last  section  of  the 

program the  printout  of  the  correspondence  of  columns  in  the  undiagonal- 

ized  and  diagonalized  matrices.  If  part  of  channel  11  is  required  for 
the  storage  of  nonzero  elements  then  certain  cells  and  commands  must 
be  altered  for  the  last  segment  of  the  program.   (Constant  in  1630, 
and  command  16^4-  n    portion).  These  could  be  changed  to  1200  wiping 
out  a  portion  of  the  program  no  longer  required.  If  there  are  less  than 
4-77  (octal)  nonzero  elements  the  program  can  be  simply  modified  to 
eliminate  the  use  of  the  tape  unit.  The  partially  ordered  P's  then  are 
stored  starting  in  0500.  The  last  iteration  is  then  not  saved. 

The  only  constants  which  must  be  filled  into  program  are  cells; 

1440  Number  of  last  cell  filled  with  a  nonzero  element 

1441  Order  of  matrix  plus  one. 
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If  P's  are  initially  completely  filled  out  start  computing  in  1200. 
If  only  the  Co  and  Ro  numbers  are  in  the  M  portion  of  the  P's  start 
computing  in  1550.  If  SEN  2020  up  program  will  type  out  all  P's  after  an 
iteration  is  completed .  SEN  204.0  must  be  up  for  program  to  recycle  after 
completing  an  iteration. 
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0000 

0000 

Print  control;  Constant  1 

1^55 

00 

0100 

0000 

0000 

Print  control 

1U56 

00 

3000 

3000 

xxxx 

Comparison  constant 

1U57 

00 

0000 

0000 

Y 

Distance  betveen  rows  (cols ) 
being  interchanged 

1U60 

00 

0000 

0000 

0030 

Shift  control  left  8  octal 

1U61 

00 

2006 

2100 

xxxx 

Comparison  constant 

1U62 

00 

X 

0000 

0000 

Variable  constant 

1U63 

00 

0000 

0000 

0006 

Shift  control  left  2  octal 

lfc6U 

00 

0000 

0000 

0100 

1U65 

10 

0100 

0000 

0000 

Print  control 

1^66 

02 

0000 

0000 

0006 

Shift  control  right  2  octal 

1U67 

00 

0000 

0000 

0010 

1U70 

00 

7700 

0000 

0000 

Extractor  for  Co 

1U71 

00 

7777 

0000 

0000 

Extractor  for  M1 

1U72 

00 

0077 

0000 

0000 

Extractor  for  Ro 

1U73 

00 

0000 

0000 

0000 

Tally  of  iterations 
I  Ir  of  last  iteration 

1U7U 

00 

0000 

0000 

0000 

1U75 

00 

0000 

0000 

0000 

£  D2  of  last  iteration 

1U76 

00 

0000 

0000 

0000 

^  D2  of  current  iteration 

1U77 

00 

0000 

0000 

0000 

X.  D2  of  current  iteration 

1500 

00 

0000 

0000 

0000 

Buffer  save  space 

1501 

00 

0000 

0000 

0000 

ti     n    tt 

1502 

00 

0000 

0000 

0000 

it     11    tt 

1503 

00 

0000 

0000 

0000 

«     it    H 

150U 

00 

0000 

0000 

0000 

it     it    it 

1505 

00 

0000 

0000 

0000 

n     H    « 

1506 

00 

0000 

0000 

0000 

n               it      it 

1507 

00 

0000 

0000 

0000 

it       11      ti 
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ADD 

INST 

M1 

M2 

M3 

REMARKS 

1510 

21 

1232 

1*55 

0001 

Carriage  return 

1511 

21 

1*73 

2100 

0003 

Print  iteration  and    D2 

1512 

21 

1232 

1*65 

0001 

Carriage  return 

1515 

17 

2020 

3000 

1525 

SEN  2020  UP,  print  P*s 

151* 

1* 

0000 

0100 

0500 

Tape  block  search  BL500 

1515 

32 

2100 

1*** 

1516 

Initialize  "buffer  cmd 

1516 

05 

3000 

3000 

xxxx 

Buffer  load 

1517 

15 

0000 

0102 

0000 

Tape  write 

1520 

35 

1516 

1*67 

1516 

Add  10  to  buffer  cmd 

1521 

3* 

1*56 

1516 

1516 

Have  we  written  enough  tape? 

1522 

17 

20*0 

3000 

1207 

SEN  20*0  UP,  recycle 

1523 

22 

0000 

0000 

0000 

Halt 

152* 

3* 

3000 

2100 

1200 

After  halt,  will  recycle  on 
compute  command 

1525 

35 

1**0 

1*52 

2001 

Form  N+l 

1526 

32 

2001 

1*** 

1527 

Extract  (N+l)  into  print  cmd 

1527 

21 

0000 

2100 

XXXX 

Print  P's 

1530 

21 

1232 

1*65 

0001 

Carriage  return 

1531 

3* 

3000 

2100 

151* 

1532 

00 

3*76 

*175 

6360 

Flexo  code 

1533 

00 

566* 

60** 

6**7 

11    tt 

153* 

00 

61*1 

7377 

66*1 

11    ii 

1535 

00 

7160 

55*5 

*731 

•1    « 

1536 

35 

'  1*76 

1*52 

1*76 

Tally  an  overflow 

1537 

35 

2100 

1*77 

1*77 

Wipe  out  overflow  bit 

15*0 

3* 

3000 

2100 

1575 

15*1 

00 

0000 

0000 

0000 

Buffer  save  space -minaccess 

15*2 

02 

0000 

0000 

0001 

Shift  control  right  1  binary 

15*3 

00 

0000 

0000 

0000 

Buffer  save  space -minaccess 

15** 

00 

1633 

0000 

0000 

Constant 

15*5 

00 

0000 

0000 

0000 

Buffer  save  space -minaccess 

15*6 

00 

0000 

0000 

0000 

15*7 

00 

0000 

0000 

0000 

Buffer  save  space-minaccess 

Only  the  Mr  portion,  i.e.,  Co  and  Ro,  need  be  initially  put  in 
the  computer.  If  only  the  n  portion  is  filled,  starting  computing 
in  cell  1550  will  complete  the  other  portions,  namely,  current  R  and 

C  and  S  for  all  elements  P.  This  portion  of  the  program  will  also 

c—  2  «sr  2 

compute  ar   initial  >_D  .  If  SEN  2010  up  this^D   will  be  entered 

and  the  initially  filled  out  P's  saved  in  tape  unit. 


ADD 

INST 

M1 

M2 

M3 

REMARKS 

COMPLETION  OF  INFORMATION 

FOR  P's:  COMPOTE  INITIALED2 

1550 

55 

2100 

2100 

1477 

Initialize  partial  ID2 
Form  N  in  M1 

1551 

30 

1440 

1U60 

1W3 

1552 

35 

1UU3 

1U5U 

14U2 

Form  N+l  in  M1 

1555 

32 

2100 

1U71 

1555 

Initialize  lookup  cmd 

1551* 

32 

2100 

14U4 

1575 

Initialize  put  away  cmd 

->1555 

35 

xxxx 

2100 

2000 

Shift  P  to  buffer 

1556 

32 

2000 

1U70 

2102 

Extract  Co 

1557 

32 

2000 

1U72 

2101 

Extract  Ro 

1560 

50 

2001 

Ikh6 

2001 

Shift  Ro  8  octal  places 

1561 

30 

2002 

1UU6 

2002 

Shift  Co  8  octal  places 

1562 

35 

2002 

2001 

2003 

>  Extract  C  and  R 

1565 

32 

2003 

1UUU 

2000 

)     into  P 

156U 

30 

2001 

IU63 

2001 

Shift  R  2  octal  places 

15*5 

35 

2002 

2001 

2001 

Form  C+R 

1566 

30 

2001 

IU63 

2002 

Shift  S  2  octal  places 

1567 

32 

2002 

1UU7 

2000 

Extract  S  into  P 

1570 

30 

2001 

1U66 

2001 

Shift  S  right  2  octal 

1571 

36 

2001 

lMn 

2001 

Form  D=S-(n+l) 

1572 

26 

2001 

2001 

2001 

Form  D2 

1573 

35 

1U77 

2001 

1U77 

Add  D2  to  partial  ^  D2 

157U 

37 

1V77 

3000 

1536 

Test  overflow 

1575 

35 

2000 

2100 

xxxx 

Put  away  P- 

1576 

35 

1555 

1U5U 

1555 

Alter  lookup 

1577 

35 

1575 

1U52 

1575 

Alter  put  away 

1— 1600 

3^ 

lW*2 

1555 

1555 

Have  we  operated  on  all  P? 
Print  ^D2 

k           l6°l 

21 

1U76 

2100 

0002 

*     1602 

32 

l440 

1U50 

2100 

Make  comparison  constant 

1603 

35 

2000 

1U6U 

2000 

for  buffer-tape  cmd 

160U 

32 

2000 

lUM* 

IU56 

I^Ofi 

17 

2010 

3000 

1U3U 

SEN  2010  UP,  shift  D2  and 
store  P»s.  DO  NOT  HAVE 

SEN  2C40  UP 

1606 

22 

0000 

0000 

0000 

CONSTANTS  OPTIMUM  SPACED 

1607 

00 

0000 

0000 

0000 

1610 

02 

0000 

0000 

0030 

Shift  control  right  8  octal 

l6ll 

00 

0000 

0000 

7700 

Extractor  for  C 

1612 

00 

0001 

0000 

0000 

Print  control;  const  1 

1613 

00 

0000 

0000 

0002 

l6lU 

00 

0000 

0000 

0077 

Extractor  for  R 

1615 

00 

0000 

0000 

7777 

Extractor  for  M3 

1616 

00 

0000 

0000 

7700 

Extractor  for  C 

I617 

00 

0000 

0000 

0006 

Shift  control  left  2  octal 

37 


FLOW  CHART  -  CORRESPONDENCE  OF  ORIGINAL  AND  NEW  COLUMNS 


j Initialize  lookup  and  put  avay 


IShift  P  to  buffer,  extract  S 
IForm  S-(n+l)    =  D  _  


\ 


Yes^ND>0 

;         \ 

\v 

(Alter  lookup 


-*No 


jPut  P  away  in  temporary  storage! 
jAlter  put  away  to  (a+l) : 


"Yes«-<  TM20  > -»No 


\ 


V 


■Have  we  tested  all  P 
i.e.,  N+l  > lookup  address 


— J" 


Initialize  Put  away  P1    to   (a-t-n) 
Form  comparison  constants 
— *f 


[initial  ize  lookup  P*,  put  ^l^y~P^] 


[Shift  P1    to  buffer  ."extract""^ 
ehACt.  P!  to  buffer,  extract  C» 


Put  P2  in 
P2  storage 

(ZZZZ 


Put  P1    in  P2  storage 
Shift  P2  to  P1  ,   S2  to  _S* 


.Alter  lookup  of  P2 
[Alter  put  away  of  P2^ 


Yes< 


Have  we  tested  all  P's 
(a+n^>lookup 


Put  P1  in  last  storage  posTt'ion '   j 
Alter  last  storage  cell  to  (a  n)-li 

,_2  ,- Is  last  storage  cell  >  (a 

> ^No - 


Form  extractor  of  C  new 

Initialize  lookup 

♦Print  heading  "Corig-Cnew" 


Shift  P  to  buffer  cell  2000 
Extract  C,  Shift  left  ?  octal 

Put  in  cell  2001 

Print  2  words-first  U  characters  only 
Alter  lookup 


Is    (a+n)>Address  of  P? 


[Add  21.2^25251525  to  itself 


No<- 


/Tesv 
-\   0'flov^>- 


— >Yes 


\. 


|  Alter  cmds  marked  * 
Print  Headings 
"Cnew-Corig" 
Print  control  to 

print  2007  and  2000 
Shifted  C  is  put  in 
I   2007 


K 

E 


j Alter  ends  marked  *  j 

back  to  original   i 

'  _sejt t i ngs  ___ \ 


jEnd  of  Program 


Figure  10 


ADD   INST 


MJ 


M2 


M3 


REMARKS 


®- 


1620 
1621 
1622 
1623 
162U 
1625 
1626 
1627 
1630 
1631 
1632 
1633 
163U 
1635 


© 


-^1636 
1637 
l6U0 
16U1 

—>l6U2 
16U3 

16U5 
16^6 

16^7 
1650 
J — 1651 
1652 
>1653 
165U 
1655 
1656 
1657 
1660 


00 
00 
00 
00 
02 
00 
00 
00 
00 
00 
00 
00 
00 
00 

30 

35 
32 
32 
35 
32 
36 
3* 
35 

35 
35 
3U 
35 
35 
36 
32 
35 
30 
32 


0000 
0000 
0000 
0000 
0000 
0000 
0000 

2525 

1100 

1736 

0000 

3*36 

6360 
6232 

lWn 

1UU7 

2100 

1630 
xxxx 

2000 

2003 

2001 
2000 

16U6 

l6te 
lW*2 
1U70 
1630 
2000 
2001 
2001 
2001 
I63O 


0000 
0777 
0000 
0000 
0000 
0000 
0000 

2525 

0000 
0000 

0777 
6232 
7327 

61*66 
1622 

2100 

1U71 
lkkk 

2100 

200U 
2007 

2100 
2100 

1U52 
1U5U 
16U2 
2100 
lMU 
1613 
lkkk 
1U52 
1U60 
lkkk 


0006 
0000 
001U 

0077 
0006 

0001 

0077 
2525 

1100 
0000 
0000 

6U77 
3136 
*551 

2007 
200U 
16U2 

l6k6 
2000 

2103 

2001 

1650- 
xxxx 

16U6 
16*2 

2002 
2000 
2001 

1675 
2001 
2003 
l*6l 


Shift  control  left  2  octal 

Extractor  for  S 

Shift  control  left  h   octal 

Extractor  for  P 

Shift  control  right  2  octal 

Extractor  for  R 

Test  constant 

Constant  for  temp  storage 

Constant 

Extractor  for  S 

Flexo  code 
11    « 

H        II 

Form  (n+l)  in  M2 

Shift  extractor  of  C  to  buffer 

Initialize  lookup 

Initialize  put  avay 

Shift  P  to  buffer 

Extract  S  into  2003 

Form  S-(n+l)«D 

Is  D  greater  than  0 

If  =0,  P  is  a  diagonal  element 

put  in  temp  storage 
Alter  put  away 
Alter  lookup 
Have  we  tested  all  P*s? 
Shift  extractor  of  Co  to  2002 
Form  1100+ (n+1) 
Form  1100-ffn-l) 
Initialize  put  away  P1 
Form  1100+n 
Shift  left  8  octal 
Make  a  comparison  constant 
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ADD   TN?T 


MJ 


M2 


M3 


FEMARKS 


1630 
2000 

16^0 
1100 

2006 
xxxx 

20C7 
2005 
20C7 
1666 
1671 
2003 
2006 
1675 
2003 
1675 
3000 
1671 
2006 
2007 
2005 
3000 
2001 
1633 
1232 
1*50 
1630 
XXXX 
2000 
2001 
2000 
171* 
1232 
2003 
1627 
1627 
1631 
1220 
1217 
3000 
2100 
15** 
1252 

1255 

0000 

3*36 
*551 

6*77 


1*5* 
1*71 
1*** 
2100 
2002 
2100 
2002 
200* 
2100 
1*5* 
1*52 
1666 
2100 
1*52 
1*5* 
l*6l 
2100 
1*** 
2100 
2100 
2100 
2100 
1*60 
1*52 
1*65 
2100 
1*71 
2100 
2002 
1*60 
1*55 
1*5* 
1*65 
171* 

1627 
3000 

1*71 

1*71 
1*** 
2100 
1627 
1*71 
1*71 
1*** 
0000 

6232 
3136 
6360 


2000 
1666 
1671 
2006 
210* 
2007 

2105 

1702-^ 

XXXX  i 

1666  ' 

1671 

1666 

XXXX 

1675 
2003  ' 

1662 

1707^-, 
1703  <r 
XXXX 

2006 

200* 

1672 

2003<- 

0003 

0001 

2002 

171* 

2000 

2101 

2001 

0002 

171* 

0001 

171* 

1627 

1731 

1710 

1717 
1716 

1653 

1627  «-* 

1710 

1717 

1716 

0000 

6h66 

6232 

7327 


Make  a  comparison  constant 
Initialize  lookup  P2 
Initialize  put  avay  P2 
Shift  P1  to  "buffer 
ExtraOt  C1  into  200* 
Shift  P2  to  buffer 
Extract  C2  into  2005 
Is  C2>C1 
Put  avay  P2 
Alter  lookup  P2 
Alter  put  avay  P2 
Tested  all  P's? 
Put  avay  P1 
Alter  put  avay  P1 
Alter  comparifon  constant 
All  P's  ordered? 

Prepare  put  avay  P1 
Put  P1  avay 
Shift  P2  to  P1 
Shift  S2  to  S1 

Prepare  a  comparison  constant 

Print  headings 

Carriage  return 

Shift  extractor  for  C  to  2002 

Initialize  lookup 

Shift  P  to  buffer 

Extract  C 

Shift  C  left  8  octal 

Print  Co  -  C 

Alter  lookup 

Carriage  return 

Have  all  C's  been  printed? 

Make  tally  constant 

Is  this  2nd  time  thru  routine? 

Change  print  control 

Change  print  control 

Change  C  put  avay  address 

Wipe  out  overflow  bit 
Change  print  control 
Change  print  control 
Change  C  put  avay  address 
Final  halt 
Plexo  code 
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